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摘要 : 传统 子 空 间 类 方法 依赖 于 阵列 相关 符 阵 ， 在 阵列 中 阵 元 出 现 失效 的 情况 下 ， 相 关 和 抵 
阵 将 会 秩 亏 ， 子 空间 类 方法 将 会 失效 。 针 对 该 问题 ， 本 文 从 协 方差 拟 合 准则 出 发 ， 将 协 方差 
矩阵 的 Toeplitz 结构 作为 约束 ， 基 于 低 秩 矩阵 重 构 的 原理 来 恢复 协 方差 矩阵 ， 而 后 采用 子 
空间 类 方法 来 进行 目标 的 方位 估计 。 在 阵 元 失效 情况 下 , 该 方法 能 有 效 重 构 阵 列 协 方差 矩阵 ， 
恢复 失效 阵 元 的 自由 度 ， 解 决 失 效 阵 元 情况 下 高 精度 目标 方位 估计 的 问题 。 数 值 仿真 表明 ， 
该 方法 在 阵 元 失效 条 件 下 , 能 够 恢复 损伤 阵列 到 正常 阵列 条 件 下 的 性 能 , 尤其 是 对 于 多 个 
标的 情况 ， 该 方法 表现 出 更 优 的 性 能 。 
关键 词 ， 阵 元 失效 ; 自由 度 恢复 ， 协 方差 矩阵 重 构 ， 低 秩 和 矩阵 台 
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High-resolution DOA estimation method based on covariance 


matrix reconstruction in presence of element failure 
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Abstract: Failed elements can cause the change of array correlation matrix. Thus, 
based on the processing of array correlation matrix, the performance of subspace 


methods would deteriorate greatly when array elements failure occurred. This paper 


presents a method to solve this problem. Starting from covariance matching criterion, 
making the Toeplitz structure of the covariance matrix as a constraint, the 
covariance matrix is recovered based on the principle of low rank matrix 
reconstruction, and then DOA estimation is conducted by the use of subspace methods. 
In case of element failure, the proposed method can reconstruct the array covariance 
matrix effectively, recover the freedom of the failed element, and solve the problem 
of high precision DOA estimation under array elements failure. Numerical simulations 
indicate that, when array element failure occurred, the method we propose have better 
performance, especially in the case of multiple targets 
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阵列 信号 处 理 是 雷达 、 声 呐 、 深 空 通 信 等 领域 中 的 关键 技术 。 在 目标 的 检测 、 
识别 与 定位 中 ， 一 般 文 献 都 假定 阵 元 是 完好 的 ， 然 而 在 实际 应 用 中 ， 由 于 恶劣 的 
自然 环境 、 人 为 干扰 以 及 阵 元 老化 、 阵 元 的 物理 损坏 等 因素 的 影响 ， 阵 元 失效 问 
题 往 往 不 可 避免 。 尤 其 在 水 下 条 件 , 环境 更 加 复杂 与 恶劣 , 如 强 水 压 、 腐 蚀 性 等 ， 
基 阵 阵 元 损坏 的 概率 更 大 ， 阵 元 失效 会 破坏 均匀 线 列 阵 的 几何 对 称 性 和 均匀 人性， 
使 阵列 输出 信号 产生 失真 ， 阵 列 响应 发 生 锐 变 ， 造 成 阵列 的 波束 图 扭曲 ， 严 重 情 
况 下 甚至 导致 阵列 失效 。 阵 元 的 失效 会 造成 检测 性 能 下 降 ， 虚 警 率 增 加 ; 在 目标 
方位 估计 中 ， 阵 元 失效 会 使 采样 协 方差 矩阵 产生 变化 ,使 得 传统 子 空间 类 的 方位 
估计 方法 失效 。 对 于 维修 复杂 、 系 统 实时 性 要 求 高 或 是 开销 昂贵 的 特殊 场合 ， 如 
何在 阵 元 失效 情况 下 恢复 阵列 的 性 能 是 一 个 需要 深入 研究 的 问题 。 

针对 阵 元 失效 问题 , 广大 学 者 分 别 从 阵 元 失效 对 于 阵列 性 能 的 影响 ， 阵 元 失 
效 所 带 来 的 旁 久 级 升 高 问题 ， 以 及 失效 阵 元 的 有 效 恢复 等 问题 展开 了 广泛 的 研 
究 。 关 于 阵 元 失效 对 于 阵列 性 能 的 影响 的 研究 , 朱德 智 等 分 析 了 均匀 线 阵 中 失效 
阵 元 在 不 同位 置 时 对 波束 主办 幅度 、 均 匀 加 窗 波束 形成 方向 图 、 锥 化 波束 形成 方 
向 图 以 及 单 脉 冲 测 角 的 影响 [1] [2]。 李 和 铮 等 分 析 均 匀 圆 阵 阵 元 失效 对 阵列 方向 图 
的 影响 , 得 出 了 阵 元 失效 的 位 置 不 同时 对 于 圆 阵 的 方向 图 的 影响 相差 不 大 的 结论 
[3] 。 针 对 和 矩形 面 阵 ， 高 飞 等 人 讨论 了 不 同 阵 元 失效 情况 对 阵列 方向 图 、 杂 波 特 
性 及 STAP 性 能 的 影响 。 从 阵列 孔径 、 杂 波 自 由 度 和 系统 自由 度 以 及 目标 相 参 积 
累 值 三 个 方面 分 析 阵 元 失效 对 空 时 自 适 应 处 理 的 影响 [4] 。 以 上 分 别 研究 了 在 均 
匀 线 阵 、 均匀 圆 阵 以 及 矩形 面 阵 下 的 阵列 性 能 影响 , 但 没有 给 出 阵 元 失效 情况 下 
的 信号 处 理 方 法 。 针 对 阵 元 失效 引起 的 波束 形成 的 旁 六 级 升 高 的 问题 , 岑 琳 等 提 
出 采用 时 间 延 迟 估计 和 线性 预测 两 种 方法 对 失效 阵 元 的 输出 进行 预测 , 抑制 旁 斩 有 
级 ， 提 高 基 阵 指向 性 [5] 。 徐 朝阳 等 提出 将 失效 阵 元 的 权重 强制 为 零 ， 根 据 期 望 
的 波束 图 对 阵列 中 剩余 正常 阵 元 的 权重 进行 重 优化 , 达到 抑制 阵列 的 旁 泊 级 的 目 
的 [6]。 针 对 阵 元 失效 的 有 效 恢复 问题 ， 杨 洋 等 以 网 格 阵列 为 例 分 析 了 阵 元 失效 
对 声场 重建 结果 的 影响 ,提出 了 忽略 缺损 阵 元 的 存在 ,把 缺损 阵列 看 作 新 的 网 格 
阵列 或 者 非 规则 阵列 ， 对 数据 进行 重组 ,利用 新 阵列 进行 重建 阵 元 失效 的 补偿 方 
法 [7] 。 张 远 彪 等 针对 基 阵 中 的 阵 元 失效 问题 ， 提 出 采用 线性 预测 方法 重 构 失 效 
阵 元 的 回 波 信号 ， 而 后 进行 合成 孔径 成 像 处 理 ， 从 而 解决 阵 元 失效 带 来 的 图 像 畸 
变 问题 [8] 。 针 对 均匀 线 列 阵 中 阵 元 失效 导致 阵列 性 能 下 降 的 问题 ， 徐 朝阳 等 利 
用 相 邻 阵 元 相对 同一 信号 源 的 接收 信号 之 间 仅 相差 一 个 固定 相 移 的 特点 , 通过 合 
成 其 他 正常 阵 元 的 输出 信号 来 恢复 失效 阵 元 的 输出 ， 抑 制 旁 为 增长 [9] 。 
Yerriswamy 等 人 [10j 在 阵 元 数据 缺损 满足 空 时 随机 分 布 的 前 提 下 ， 提 出 通过 低 
秩 和 矩阵 填充 的 方法 来 恢复 阵列 在 空 时 上 的 缺损 数据 , 然后 在 数据 域 直接 采用 矩阵 
束 的 方法 进行 目标 方位 估计 。 在 实际 中 , 假设 的 空 时 数据 缺损 的 随机 分 布 条 件 往 
往 不 能 同时 满足 ， 阵 元 损坏 会 导致 某 一 空间 位 置 无 法 接收 到 任何 数据 ， 从 而 使 得 
该 方法 失效 。 杨 东 等 人 [11] 考虑 了 空间 采样 中 某 一 阵列 完全 损坏 或 者 未 采样 情 
况 下 的 问题 ， 提 出 利用 等 距 线 阵 数据 间 的 等 比特 性 ,通过 将 某 一 时 间 采 样 上 的 数 
据 重 排 ， 使 得 数据 满足 随机 分 布 条 件 ， 进 而 将 问题 等 效 为 一 低 秩 矩 阵 填 充 问题 。 
该 方法 虽然 解决 了 某 一 阵 元 损坏 或 者 降 采 样 时 的 空 时 数据 非 随机 分 布 ,但 存在 2 
个 重要 缺陷 ， 首 先 ， 数 据 重 排 后 的 矩阵 填充 是 个 降 维 处 理 ， 又 要 保证 低 秩 条 件 ， 
极 大 地 降低 了 阵列 的 自由 度 。 其 次 , 逐个 快 拍 进行 数据 填充 所 需 的 计算 量 非 常 大 ， 
在 实际 应 用 中 不 太 现 实 。 在 阵 元 失效 的 情况 下 如 何 有 效 恢复 阵列 信号 ,进行 高 分 
辨 的 方位 估计 、 成 像 [12] 已 成 为 阵列 信号 处 理 中 的 一 个 重要 研究 内 容 


针对 均匀 线 列 阵 的 阵 元 失效 问题 ,在 假设 阵 元 损伤 为 随机 分 布 ， 阵 元 失效 个 
数 不 超 过 一 半 的 情况 下 , 本 文 从 协 方差 拟 合 准则 出 发 , 以 协 方差 矩阵 的 Toeplitz 
结构 作为 约束 条 件 [13], 利用 协 方差 矩阵 的 低 秩 特性 来 恢复 协 方差 矩阵 ， 而 后 采 
用 子 空间 类 方法 进行 目标 的 方位 估计 。 在 阵 元 失效 情况 下 ， 该 方法 能 有 效 恢 复 协 
方差 矩阵 ,恢复 失 效 阵 元 的 自由 度 ,， 解决 失效 阵 元 情况 下 高 精度 目标 方位 估计 问 
题 。 
2 均匀 线 列 阵 模型 
2.1 无 损伤 均匀 线 列 阵 楼 型 

考虑 包含 有 M 个 阵 元 的 均匀 线性 阵列 。 假 设 有 兢 个 非 相 关 窗 带 目标 信 
号 分 别 从 远 场 gp,py…prx 方 位 照射 到 阵列 。 不 同时 间 快 拍 表 示 为 s.00)， 
t=1,2,…,T，k=L2,…,K 。 那 么 该 阵列 接收 的 信号 矢量 可 以 表示 为 : 


x()= Da (pd) + nd) = As(O) + nd) (1) 
k=l 


TS 


j25 sin' j2 sin(@, F ~、 攻 
式 中 : a(p)=| he2 ,e249| ec™m 为 gp 对 应 的 导向 矢量 ，d 表 示 阵 


元 位 置 ， ad=dG1) ，i=l2M 。 对 pw 进行 归 一 化 ， 令 


A=[a(@),a(8),…,a(9gx)| 为 阵列 流 型 , 假定 sG 服从 方差 为 的 独立 高 斯 分 
布 ， 表 示 为 sb =[s(D,s(，sl 。 噪 声 n( 力 为 服从 复 高 斯 分 布 NC(0,0?) 的 
独立 同 分 布 随 机 变量 ,均值 为 0, 方差 为 o? 。 那 么 信号 x( 人 的 协 方差 算 
阵 可 以 表示 为 : 


of 
R, = E[x(D)x" (0)]=A 9 . A +ol=ARA" +ol=T()+diag(o’) (2) 
Ok 
式 中 : R=E[s(0)s*(t)]=diag([o?,.…o2]))eR** 为 信 源 的 协 方差 矩阵 ; ec? 为 第 K 个 信 
源 的 功率 ; I 表示 单位 矩阵。 

6 -> 09， k=1,2,..,K (3) 
实际 中 ， 协 方差 矩阵 是 未 知 的 ， 只 有 在 采样 时 间 无 限 长 的 情况 下 才能 

够 精确 地 估计 。 一 般 协 方差 矩阵 可 以 通过 如 下 的 采样 协 方差 矩阵 来 获得 : 
良 _ -了 XO) (4) 


在 阵列 无 损伤 情况 下 ， 可 以 采用 子 空间 类 方法 如 MUSIC、ESPRIT 等 来 
通过 寻找 伪 谱 峰值 所 对 应 的 角度 来 进行 目标 方位 估计 ， 但 是 在 阵 元 存在 损 
伤 的 情况 下 ， 会 导致 秩 亏 以 及 相位 模糊 ， 使 得 传统 子 空间 类 方法 失效 。 本 文 
从 协 方差 拟 合 准则 出 发 ， 以 协 方差 矩阵 的 Toeplitz 结构 作为 约束 条 件 ， 利 用 协 
方差 矩阵 的 低 秩 特征 来 恢复 协 方差 矩阵 , 而 后 采用 子 空 间 类 方法 来 进行 目标 的 方 
位 估计 。 

2.2 损伤 均匀 线 列 阵 模 型 

这 里 我 们 作 如 下 假定 : 

(1) 已 知 损伤 阵 元 的 空间 位 置 。 

(2) 这 里 所 指 的 阵 元 失效 是 指 阵 元 完全 故障 ， 接 收 不 到 任何 幅度 和 相 


位 的 信息 ， 也 就 是 说 失效 阵 元 没有 可 以 利用 的 任何 信息 。 若 是 第 7 个 阵 元 
失效 ， 那 么 导向 矢量 a(g) 中 的 e234” 将 被 0 所 替代 。 


(3) 损伤 阵 元 的 空间 位 置 是 随机 分 布 的 ， 且 阵 元 损伤 个 数 不 超 过 一 半 。 
损伤 均匀 线 列 可 以 看 做 是 无 损伤 均匀 线 列 阵 的 空间 域 稀疏 采样 ， 也 就 
说 损伤 阵列 输出 可 以 看 做 是 原始 均匀 线 阵 输出 的 一 个 子 集合 。 假 若 损伤 阵 


元 的 阵 元 位 置 集合 为 YY， 那么 VcIM]， 我 们 认为 该 集合 指标 按照 升序 排列 ， 


= 为 除 损伤 阵 元 的 阵列 阵 元 的 个 数 ， 那 么 损伤 阵列 的 导向 矢量 可 以 表示 
为 


.2 & .27T i: 
jd sin(@; ) j——dy, sin(9, ) 


oy(Gi)=les et eC™ (5) 


假定 损伤 阵 元 位 置 已 知 ， 定 义 一 个 选择 矩阵 Ty ef0,1}”， 其 中 工 , 的 第 j 行 除 


了 第 个 元 素 为 1， 其 余 全 部 为 零 。 那 么 损伤 阵 元 的 导向 矢量 可 以 表示 为 


au (0.)=Ta(9) (6) 
阵列 输出 可 以 表示 为 
xm(D = 4os(D +nao(D =To4s(D +Ton(®) (7) 
那么 此 时 的 协 方差 矩阵 为 
Rs, = lim xn(Dxg(D =TyT(WTY + Tydiag(o TY ST(w) +diag(o?) (8) 


这 里 Tj,(w)aTyT(WrT5 ， diag(o2) 和 Twydiag(o3T% 。 实 际 中 ， 采 用 采样 协 方差 表示 为 
入 四 1 y 家 
及 i Xim (Dxin (t) o 


3 ”基于 低 秩 矩阵 约束 的 协 方差 矩阵 重 构 
在 无 损伤 阵列 模型 下 ， 令 B=A)Rs4e(o) ， 很 容易 得 出 B>0， 并 且 


rank(B)=K <M-!1l,， 进一步 


K K 
Bi = > caic(w Jo (9)= Yee (9) 
k=1 k=1 
B 是 一 个 (Hermitian) Toeplitz 和 矩阵， 可 表示 为 B=T(u)，ueC*， 且 有 
ul L2 Um 
的 | > (10) 
U Um ul 


在 采样 协 方差 矩阵 两 边 减 去 协 方差 矩阵 ， 可 以 得 到 


E=R, -R = 入 -TO)-diag(o2) (11) 


式 中 ，E 是 测 量 误差， 因为 快 拍 数 的 有 限 ， 使 得 该 测量 误差 不 可 能 为 零 ， 其 
中 包含 本 信 筷 与 信号 ， 信 和 号 与 噪声 之 间 的 交叉 相关 项 ， 随 着 快 拍 数 的 增加 ， 
该 误差 会 逐渐 变 小 。 

令 D={m-mt+tl:im,m et,m2>m}c[M]; 当 了 定义 在 均匀 线 列 阵 上 的 ，D 为 元 
余 阵 ， 且 损伤 阵列 可 视 为 非 均 匀 线 阵 上 的 合作 阵 。 为 保证 损伤 阵列 的 元 余 性 ， 假 
定 损伤 阵 元 位 置 是 随机 的 且 个 数 不 超过 一 半 ,， 这 种 假设 在 实际 中 是 合理 的 。 该 假 
设 是 充分 非 必要 条 件 ， 即 使 合作 阵 不 是 元 余 阵 ， 也 可 以 保证 协 方差 的 重 构 ” 。 

对 于 均匀 线 阵 我 们 可 以 通过 低 秩 矩阵 重 构 理 论 ， 通 过 以 下 方式 来 重 构 
协 方差 矩阵 


[TQ)] st Eb <p,T(u)>0 (12) 


这 里 ，o” > 0 意味 着 非 负 向 量 。p 是 一 个 用 户 自 定义 的 上 界 ， 这 里 采用 协 方差 拟 
合 的 策略 来 确定 该 界限 ， 首 先 通 过 向 量化 E 满 足 渐进 正 态 分 布 
vec(E) ~ ASN(0, W) 其 中 ，W =1/N(R'@R)，@ 为 Kronecker 矩阵 积 。 实 际 中 ， 
可 以 通过 估计 得 到 W =1/N(R @R)， 那 么 可 以 得 到 


W vec(E)~ AsN(0,1,,) (13) 
直接 的 结果 是 
ve ~ Asz2(OV2) (14) 
Asx (N ) 是 自由 度 为 N” 的 卡 方 分 布 ， 引 入 一 个 参数 7， [ 必要 那么 
通过 如 下 的 秩 最 小 问题 来 求解 协 方差 矩阵 重 构 的 问题 
| St W 3vec(E) | <n,T(u)>0 (15) 


但 是 ， 该 问题 是 一 个 NP 难 问 题 ， 因 此 采用 凸 松弛 将 伪 秩 范 数 松弛 到 原子 核 范 数 
或 者 正 半 定 和 矩阵 的 迹 范 数 。 因 此 该 优化 问题 可 以 表示 为 


2 


eal St Wzvec(E) <n,T(u)>0 (16) 


一 旦 将 优化 解 立 得 到 ， 那 么 通过 T(u) =T(O) 得 到 重 构 协 方差 矩阵 ， 进 一 步 可 
以 通过 root-MUSIC 估计 得 到 目标 方位 。 因 为 T(0) 上 共有 Toeplitz 结构 ， 也 可 以 
通过 范 德 蒙特 分 解 来 有 效 估计 目标 方位 。 

当 阵 列 存在 损伤 时 ， 阵 列 协 方差 重 构 问题 可 以 描述 为 


2 


Wy2vec(E,)| ~ Asx’(L) (17) 
采用 如 下 的 迹 最 小 问题 来 重 构 协 方差 
| Ss.t IW,?vec(E,) <7,T(u)>0 (18) 


7 是 自由 度 为 损伤 的 阵 元 个 数 上 的 渐进 卡 方 分 布 。 该 方法 能 够 潜在 的 恢复 出 来 
损伤 阵 元 的 协 方 差 矩阵 。 因 此 整个 算法 的 流程 为 ; 


协 方差 矩阵 重 构 目 标 方位 估计 算法 流程 : 
计算 损伤 阵列 的 协 方差 矩阵 R ; 
: 计算 加 权 权 重 W ， 以 及 协 方差 拟 合 误差 参数 7 ; 
: 求解 迹 最 小 化 问题 〈18)， 重 构 协 方差 矩阵 T(u) ; 
4: 采用 求 根 MUSIC 求解 损伤 阵 元 的 目标 方位 
4 仿真 与 分 析 
本 节 对 基于 协 方差 矩阵 重 构 的 高 分 辨 佑 计 方法 进行 仿真 分 析 ， 以 验证 
算法 的 有 效 性 。 在 下 面 所 有 仿真 中 均匀 线 阵 阵 元 个 数 均 为 10 个 ， 阵 元 间隔 
为 半 波 长 ， 目 标 为 远 场 非 相干 窄带 目标 。 噪 声 为 均值 为 0， 方 关 为 1 的 高 斯 
白 噪声 。 其 中 阵列 中 有 4 个 阵 元 以 随机 分 布 的 形式 损坏 ， 在 均匀 线 列 阵 中 分 
别 出 现 的 位 置 为 [2，4，6，9]j。 破 损 的 阵 元 位 置 如 图 1 所 示 ， 其 中 〇 表示 损 
伤 阵 元 指标 ，X 表示 损伤 阵 元 指标 ， 满 足 阵 元 损伤 数目 不 超过 阵 元 数目 的 
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Fig.1 Array Configuation 


实验 1: 阵 元 损伤 情况 下 的 性 能 比较 


ar 


中 


全 


[a 


Normalized power spectrum 
已 [= [= © [= 
a 


一 一 元 好 阵列 MUSIC 


阵 页 
一 一 一 损伤 阵列 CMR Root MUSIC 
…… 真实 DOA 


0 ， 
-20 -15 -10 -5 5 10 15 20 
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图 2 不 同情 况 下 的 空间 谱 比 较 
Fig.2 Spatial Spectrum of all methods 
图 2 是 在 两 个 目标 方位 分 别 为 [-4”, 4”]。 可 以 看 出 , 将 损伤 阵 元 输出 
看 作为 0 时 ， 采 用 MUSIC 方法 无 法 分 辨 出 两 个 近 空 间 目 标 ， 但 是 通过 重 构 
后 协 方差 采用 求 根 MUSIC 可 以 达到 与 完好 阵列 接近 的 性 能 。 说 明 该 方法 在 
恢复 阵列 自由 度 上 的 有 效 性 。 
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实验 2: 测 向 偏差 与 角度 分 隔 的 关系 

实验 条 件 : 输入 信 噪 比 为 10dB。 假 设 有 2 个 等 功率 的 独立 非 相 干 信号 
同时 入 射 到 10 元 均匀 线 列 阵 中 ， 进 行 200 次 蒙特 卡 洛 实验 ， 分 析 2 个 目标 
的 平均 测 向 偏差 。 一 个 目标 的 角度 固定 为 -10” ， 另 一 个 目标 的 角度 间隔 从 
以 0.5” 步 长 远离 该 目标 ， 信 噪 比 固定 为 0dB， 快 拍 数 为 500 次 ， 分 别 在 阵 
元 无 损伤 采用 MUSIC 方法 ， 阵 元 损伤 情况 下 采用 MUSIC 方法 以 及 阵 元 损伤 
情况 下 采用 协 方差 矩阵 重 构 的 方法 来 重新 获得 协 方差 矩阵 , 采用 求 根 MUSIC 
方法 来 估计 方位 ，MUSIC 方法 中 对 [-90”, 90”] 空 域 中 以 1 为 间隔 进行 扫 
描 ， 验 证 该 方法 的 测 向 偏差 与 角度 分 隔 的 关系 ， 如 图 3 所 示 。 

在 随机 出 现 4 个 阵 元 损伤 的 情况 下 ， 采 用 协 方差 矩阵 重 构 结合 求 根 
MUSIC 方法 能 够 在 目标 角度 很 小 的 情况 下 得 到 很 好 的 估计 性 能 , 甚至 优 于 完 
整 阵列 情况 下 MUSIC 的 估计 性 能 ， 其 原因 是 一 方面 采用 Toeplitz 结构 来 约 
束 协 方差 窍 阵 ， 另 一 方面 是 因为 采用 了 求 根 MUSIC， 人 性 能 比 MUSIC 要 更 优 。 
在 阵列 损伤 情况 下 ，MUSIC 方法 在 角度 分 隅 较 远 时 ， 也 可 以 估计 出 目标 ， 因 
为 损伤 阵列 中 ， 依 然 存 在 两 个 半 波 长 的 阵 元 ， 但 是 对 于 多 个 目标 ， 其 性 能 
就 会 显现 出 来 。 采 用 协 方差 矩阵 重 构 结合 求 根 MUSIC 方法 在 两 个 目标 间隔 
为 6"” 时 ， 估 计 偏 差 已 经 逐渐 接近 为 零 ， 而 完好 阵列 在 两 个 目标 间隔 为 8° 
时 ,估计 偏差 才 副 近 为 0， 因 此 该 方法 可 以 很 好 地 恢复 出 完好 阵列 的 估计 性 
能 ， 这 正 是 该 方法 的 目的 所 在 。 
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图 3 估计 偏差 与 角度 分 离 的 关系 曲线 


Fig.3 Bias of two source separation for the different case 


实验 3: 验证 算法 在 多 目标 情况 下 的 DOA 估计 性 能 
采用 均 方 根 误差 (Root Mean Square Error，RMSE) 作为 性 能 衡量 准 
则 ， 定义 为 : 


1 
RMSE = | 一 gm) (14) 


式 中 ; M 为 给 定 场景 下 的 蒙特 卡 洛 仿真 次 数 ; k 为 入 射 信号 个 数 ， 名 和 分 
别 为 第 m 次 试验 中 第 个 信号 方向 的 真实 值 和 估计 值 。 
假设 空间 远 场 非 相干 窄带 目标 个 数 分 别 为 5，6，7，8 分 别 均匀 分 布 在 
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[-60" ,60" ]， 图 4、 图 5 分 别 给 出 了 RMSE 与 快 拍 数 、 信 噪 比 的 关系 变化 
曲线 。 图 4 是 在 信 噪 比 为 0dB 条 件 下 ， 蒙 特 卡 洛 次 数 为 200 次 ， 多 目标 的 
均 方 根 误差 与 快 拍 数 关 系 曲线 。 可 以 看 出 在 随 着 目标 个 数 的 增加 ， 估 计 精 
度 相 对 的 降低 ， 但 是 在 目标 数 为 5，6，7 的 情况 下 ， 该 方法 仍 能 够 具有 较 
高 的 精度 。 


RMSE(") 


Snapshot 
图 4 RMSE 与 快 拍 数 的 关系 变化 曲线 


Fig.4 The RMSE against different Snapshot of multiple source for different methods 


图 5 是 在 快 拍 数 为 500 次 的 条 件 ， 蒙 特 卡 洛 次 数 为 200 次 ， 多 目标 的 均 方 
根 误 差 与 信 噪 比 关 系 曲线 。 可 以 看 出 在 小 于 0dB， 该 方法 精度 较 低 ， 但 是 随 
着 信 噪 比 的 增加 ， 目 标 数 分 别 为 5，6，7 的 情况 下 ， 仍 然 具 有 较 好 的 分 辩 精 
度 ， 因 为 目标 损伤 的 个 数 为 4 个 ， 传 统 的 MUSIC 方 法 因为 存在 阵 元 缺损 ， 不 
能 分 辨 超过 阵 元 个 数 的 目标 
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图 5 RMSE 与 信 噪 比 的 关系 变化 曲线 
Fig.5 The RMSE against different SNR of multiple source for different methods 
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